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Abstract 

Objective: Several studies have reported the association between methylenetetrahydrofolate reductase (MTHFR) C677Tand 
A1298C polymorphisms and nonsyndromic cleft lip with or without palate (NSCL/P) in Asian populations. However, findings 
have been conflicting. In order to investigate the association, a meta-analysis was performed. 

Methods: We searched Pubmed, MedLine and EmBase database to selected eligible studies. The pooled odds ratios (ORs) 
with 95% confidence intervals (95%Cls) were calculated using fixed effects model or random effects model to assess the 
association between MTHFR polymorphisms and NSCL/P in both Asian children and mothers. 

Results: Finally, nine case-control studies were included. Overall, the MTHFR C677T polymorphism and NSCL/P showed 
pooled ORs (95%CI) of 1.41(1.23-1.61) in Asian children, and 1.70(1.19-2.42) in Asian mothers. Subgroup analyses by 
geographical locations further identified the association in Eastern Asian children, Western/Central Asian children and 
mothers, but not in Eastern Asian mothers. However, no significant relationship between MTHFR A1298C polymorphism 
and NSCL/P was found in this meta-analysis. 

Conclusions: The MTHFR 677T allele was associated with an increased risk of NSCL/P in Asian populations. 
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Introduction 

Nonsymdromic cleft lip with or without palate (NSCL/P) is one 
of the most common congenital malformations of the head and neck 
area in the world, occurring in approximately 1 in every 700 live 
births [1]. Around the globe, the distribution of prevalence rate of 
NSCL/P is not homogeneous, and higher values are found in parts 
of Asia (China, Japan) and Latin America, whereas Israel, South 
Africa and Southern Europe showed the lowest values [2]. The 
etiology of NSCL/P is considered to be multifactorial, causing by 
both hereditary factors and environmental factors. Epidemiologic 
studies revealed that mothers who used multivitamins containing 
folic acid showed a lower risk of having an offspring with NSCL/P 
compared with mothers did not use multivitamins [3]. However, it is 
still unknown which ingredients in multivitamins contribute to this 
risk reduction. It is widely accepted that folic acid is a key factor in 
the development of craniofacial structures. Several studies have 
reported a reduced risk of NSCL/P when mothers used either folic 
acid supplements or dietary folate during pregnancy [4,5]. Some 
other studies, however, have provided variable or even contradic- 
tory results [6-8] . Therefore, genetic variations in folate metabolism 
gene were believed to affect individual susceptibility to NSCL/P. 



Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme 
in the metabolism of folate. This enzyme irreversibly catalyzes the 
reaction of 5, 1 0-methylenetetrahydrofolate to 5-methyltetrahy- 
drofolate, which is the primary circulating form of folate [6] . The 
most common single nucleotide polymorphism (SNP) is MTHFR 
C677T, which results in an alanine to valine exchanges and is 
associated with reduced enzyme activity. MTHFR A1298C, 
another common functional polymorphism, also affects enzyme 
activity [7,8]. The two SNPs have been studied as candidate 
genetic factors for NSCL/P risk [9-12]. 

The first report evaluating the role of MTHFR gene 
polymorphisms in the development of NSCL/P was conducted 
by Tolarova et al.(1998) [13]. The findings suggested that 
MTHFR C677T variant genotype increased the risk of NSCL/ 
P. They also found a significant association between MTHFR 
A1298C polymorphism and NSCL/P. Since then, a great number 
of studies have been conducted, but the results were inconsistent 
[14-16]. In order to elucidate the role of MTHFR C677T and 
A1298C polymorphisms in NSCL/P, several meta-analyses were 
performed [17-20], but the relationships especially among Asian 
subjects remained unclear. Recently, some new case-control 
studies with large sample size have been published [21-23]. 
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61 potentially related articles identified 
from PubMed, MedLine and EmBase 
database (duplications excluded) 



Full-text reviewed 



Abstract screened 



22 articles excluded: 

Not for MTHFR C677T or 

A1298C(n=8), 

Review or Meta-analysis or Letter 
articles (n=14) 



3 1 excluded: 

Not case- control design (n=4) 
Not Asian populations (n=27) 



9 studies from 8 articles included 



8 studies for C677T: 
8 for children 
6 for mothers 



6 studies for A1298C: 
6 for children 
3 for mothers 



Figure 1. Flow chart explaining the selection of eligible studies included in the meta-analysis. 

doi:1 0.1 371 /journal.pone.0088242.g001 



Hence, we performed this meta-analysis to derive a more precise 
estimation of the association of MTHFR C677T and A1298C 
polymorphisms with NSCL/P in Asian populations. 

Materials and Methods 

Literature search and selection 

Eligible literatures were screened from PubMed, MedLine and 
EmBase database (up to May 31, 2013). We used the following 
keywords and subject terms: "methylenetetrahydrofolate reduc- 
tase" or "MTHFR" and "cleft lip" or "cleft palate". Furthermore, 
we manually searched references in the eligible articles. There was 
no language limitation. And the search results were limited to 
humans. 

Criteria of inclusion and exclusion 

Inclusion criteria were showed as following: (1) evaluation the 
association between MTHFR C677T and/or A1298C polymor- 
phisms and NSCL/P risk in Asian populations; (2) case-control 
studies; (3) the frequencies of all genotype distribution or other 
available data for estimating the OR (95% CI). If overlapping 
cases or controls were presented in multiple studies, the most 
recent publication or the largest study was included. The exclusion 
criteria were: (1) not using a case-control study design; (2) had no 
detailed data on genotype distribution; (3) not in Asian popula- 
tions. 

Data extraction 

The following information was carefully and independendy 
collected from each eligible study by two authors: first author's 
name, year of publication, country, population of cases and 
controls, genotyping method, and numbers of each genotype. The 
difference was settled by reaching an agreement among all 
authors. 



Statistical analysis 

To assess the study quality, Hardy-Weinberg equilibrium 
(HWE) in controls in each study was calculated by chi-squared 
test. P value<0.05 was considered a departure from HWE. The 
associations between MTHFR C677T and A1298C gene poly- 
morphisms and NSCL/P risk were estimated by the odds ratios 
(ORs), together with the 95% confidence interval (95%CI). The 
significance of the pooled OR was determined by the Z test, with 
P<0.05 considered significant. Heterogeneity between studies was 
assessed by 0_ test. If P<0.05, the heterogeneity was considered 
statistically significant. The I 2 values were used to quantify the 
percentage of the total variation among studies when heteroge- 
neity was assessed. The I value ranged from 0 to 100%. 25%, 
50%, 75% expressed low, moderate, high heterogeneity, respec- 
tively. When I 2 <50%, a fixed effects model was applied to 
estimate the pooled results. Otherwise, the random-effect model 
was used. 

Publication bias was investigated visually in a funnel plot of log 
(OR) against its standard error (SE). An asymmetric plot suggested 
possible publication bias. And the degree of asymmetry was 
assessed by Egger's test. Sensitivity analysis was performed by 
omitting each study in turn to assess the results stability. 

Software ST ATA version 1 1.0 was used for all analyses. 

Results 

Selection of studies 

A flow chart summarizing the process of study selection is 
shown in Figure 1. Based on the inclusion and exclusion criteria, 
eight articles were included after full-text reviewing [14,16,21-26]. 
Worthy of note, a study by Aida et al.(2012) consisted of two 
groups of control was treated as two studies [21]. Thus, a total of 
nine studies from eight articles were included in this meta-analysis. 
Among them, nine studies provided date on MTHFR C677T 
polymorphism (eight studies were carried out in Asian children 
and six in Asian mothers), and 6 studies on A1298C polymorphism 
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Study 
ID 



Events. Events, % 
OR (95% CI) cases controls Weight 



Kumah et al. (2013) 

Aida et al. (2012) 

Aida et al. (2012) 

Han et al. (2011) 

Guo et al. (2009) 

AN et al. (2009) 

Wan et al. (2006) 

Shotelersuk et al. (2003) 

Overall (l-squared = 42.3%, p = 0.096) 



.274 



1.50 (1.15,1.95) 155/934 110/938 25.01 

1.22 (0.71,2.10) 34/112 40/152 6.46 

1.09 (0.65.1.83) 34/112 53/186 7.58 

1.37 (1.03. 1.81) 176/374 168/426 22.72 

1.07 (0.72.1.58) 101/192 105/206 13.12 

1.97 (1.34.2.89) 109/646 40/428 10.93 

■ 2.22 (1.35. 3.65) 77/152 38/120 5.72 

0.92 (0.55.1.53)25/218 50/404 8.48 

1.41 (1.23,1.61) 711/2740 604/2860 100.00 



3.65 



Figure 2. Forest plot of association between MTHFR C677T polymorphism and NSCL/P risk in children (T vs. C). 

doi:1 0.1 371 /journal. pone.0088242.g002 



Study 
ID 



Events. Events. % 



OR (95% CI) cases controls Weight 



Aida et al. (2012) 

Aida et al. (2012) 

Wang et al. (2012) 

Guo et al. (2009) 

Ali et al. (2009) 

Shotelersuk et al. (2003) 

Overall (l-squared = 69.5%, p = 0.006) 



NOTE: Weights are from random effects analysis 



.255 




2.33 (1.38,3.92) 49/108 40/152 15 77 

2.08 (1.27,3.42) 49/108 53/186 16 37 

2.24 (1.41,3.57) 117/178 59/128 17.00 

0.83 (0.56,1.24) 93/194 107/204 18.57 

1.96 (1.22,3.15) 39/232 40/428 16.82 

1.47 (0.86.2.51) 23/134 50/404 15.47 

1.70 (1.19,2.42) 370/954 349/1502 100.00 



3.92 



Figure 3. Forest plot of association between MTHFR A1298C polymorphism and NSCL/P risk in mothers (T vs. C). 

doi:1 0.1 371 /journal.pone.0088242.g003 
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Table 2. Overall and subgroup results of the association between MTHFR C677T polymorphism and NSCL/P. 





Study(n) 


T vs. C 


TT vs. CC 




CT vs. CC 


CT/TT vs. CC 


TT vs. CC/CT 




OR(95%CI) P 


OR(95%CI) 


P 


OR(95%CI) P 


OR(95%CI) P 


OR(95%CI) P 


Overall 


Children(n = 8) 


1.41(1.23-1.61) <0.001 


1.85(1.30-2.63) 


0.001 


1.61(1.21-2.14) 0.001 


1.61(1.22-2.11) 0.001 


1.36(1.00-1.86) 0.052 


Mothers(n = 6) 


1.70(1.19-2.42) 0.003 


2.58(1.09-6.09) 


0.031 


1.66(1.08-2.55) 0.021 


1.83(1.17-2.88) 0.008 


1.95(1.32-2.86) 0.001 


Eastern Asia 


Children(n = 4) 


1.31(0.95-1.80) 0.099 


1.79(1.16-2.76) 


0.009 


1.68(0.88-3.22) 0.117 


1.66(0.89-3.10) 0.111 


1.26(0.87-1.83) 0.217 


Mothers(n = 3) 


1.39(0.75-2.54) 0.293 


2.19(0.46-10.44) 


0.326 


1.09(0.72-1.65) 0.688 


1 .28(0.64-2.59) 0.487 


2.06(0.64-6.61) 0.224 


Western/Central 
Asia 


Children(n = 4) 


1.50(1.25-1.81) <0.001 


1.96(1.08-3.55) 


0.026 


1.56(1.25-1.95) <0.001 


1.61(1.30-1.99) <0.001 


1.44(0.53-3.95) 0.477 


Mothers(n = 3) 


2.11(1.58-2.81) <0.001 


3.66(1.78-7.52) 


<0.001 


2.40(1.64-3.50) <0.001 


2.53(1.75-3.64) <0.001 


2.38(1.22-4.66) <0.001 
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(six in Asian children and three in Asian mothers). The main 
characteristics of selected studies were summarized in Table 1 . 

Meta analysis 

The association between MTHFR C677T and NSCL/P risk 
was shown in Table 2. Significant heterogeneity between studies 
was found in the majority of the genetic comparisons, except for 
allele comparison (T vs. C: I 2 = 42.3%), homozygote comparison 
(TT vs. CC: I 2 = 28.3%) and recessive comparison (TT vs. CC/ 
CT: I = 17.4%) in Asian children, and recessive comparison (TT 
vs. CC/CT: I 2 = 49.8%) in Asian mothers. OveraU, the variant T 
allele of MTHFR C677T increased risk of NSCL/P, when 
compared with the wild-type C allele. The OR (95% CI) was 
1.41(1.23-1.61, P<0.001, Figure 2) in Asian children, and 
1.70(1.19-2.42, P = 0.003, Figure 3) in Asian mothers, respective- 
ly. In the subgroup analyses of geographical location, the similar 
associations were also observed in Eastern Asian children, 
Western/ Central Asian children and mothers, but not in Eastern 
Asian mothers. 

For MTHFR A1298C polymorphism, the main results were 
shown in Table 3. Sound homogeneity between studies was seen, 
except for allele comparison (C vs. A: I 2 = 71.0%), heterozygote 
comparison (AC vs. AA: I = 76.0%) and dominant comparison 



significant associations between MTHFR A1298C polymorphism 
and NSCL/P were found under any of genetic model, even in the 
subgroup analyses by geographical location, neither in Asian 
children nor in Asian mothers. 

Sensitivity analysis and publication bias 

Elimination of each study made no qualitative difference on the 
pooled OR values, which indicated that the final results of the 
meta-analysis were stable. The publication bias of the included 
studies was assessed by the Funnel plot and Egger's test. The 
shapes of Funnel plot in each genotype comparison indicated no 
obvious asymmetry (Figure 4). And Egger's test provided statistical 
evidence for the funnel plot symmetry. No significant publication 
bias was found in the studies. 

Discussion 

In this meta-analysis, we investigated the association between 
MTHFR C677T and A1298C polymorphisms and NSCL/P, with 
a total of nine case-control studies included. The pooled results 
indicated that there was an obvious association between MTHFR 
C677T polymorphism and NSCL/P risk in Asian children under 
all models: allele contrast (T vs. C), homozygote (TT vs. CC), 
heterozygote (CT vs. CC) and dominant (TT/CT vs. CC) models, 



(AC/CC vs. AA: V = 74.8%) in Asian children. However, no 
Table 3. Overall and subgroup results of the association between MTHFR A1298C polymorphism and NSCL/P. 





Study(n) 


T vs. C 




TT vs. CC 




CT vs. CC 




CT/TT vs. CC 




TT vs. CC/CT 






OR(95%CI) 


P 


OR(95%CI) 


P 


OR(95%CI) 


P 


OR(95%CI) 


P 


OR(95%CI) 


P 


Overall 


Children(n = 6) 


1.03(0.79-1.35) 


0.819 


1.08(0.78-1.48) 


0.652 


1.01(0.68-1.48) 


0.982 


1.03(0.71-1.48) 


0.881 


1.09(0.80-1.48) 


0.578 


Mothers(n = 3) 


0.94(0.73-1.20) 


0.593 


0.78(0.40-1.50) 


0.452 


0.98(0.71-1.35) 


0.905 


0.95(0.70-1.30) 


0.751 


0.79(0.41-1.50) 


0.467 


Eastern Asia 


Children(n = 3) 


1.05(0.61-1.82) 


0.855 


0.97(0.47-1.99) 


0.930 


1.16(0.44-3.06) 


0.767 


1.13(0.49-2.61) 


0.768 


0.96(0.48-1.93) 


0.907 


Mothers(n = 1) 




















0.224 


Western/Central Asia 


Children(n = 3) 


1.04(0.73-1.49) 


0.819 


1.23(0.53-2.85) 


0.626 


0.94(0.76-1.16) 


0.551 


0.98 (0.66-1.44) 


0.901 


1.23(0.60-2.53) 


0.568 


Mothers(n = 2) 


0.83(0.62-1.12) 


0.218 


0.81(0.20-3.27) 


0.767 


0.81(0.54-1.20) 


0.283 


0.79(0.54-1.15) 


0.218 


0.88(0.24-3.17) 


0.839 
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/ s 

/ • \ 




/ • V 

1 1 1 r 

-2 0 2 4 

OR 



Figure 4. Funnel plot analysis to detect publication bias. 
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except for recessive (TT vs. CT+CC) model. Meanwhile, the 
variant T allele of Asian mothers was significantly associated with 
a 1.70-fold increased risk of having a NSCL/P offspring. 
Subgroup analyses by geographical location further identified this 
association in Eastern Asian children, Western/Central Asian 
children and mothers, but not in Eastern Asian mothers. 
Regarding MTHFR A1298C polymorphism, the pooled results 
revealed that A1298C was not related to children's or mothers' 
NSCL/P susceptibility under any of the genetic model, even in the 
subgroup analyses by geographical location. The above results 
suggested that MTHFR C677T polymorphism was a risk factor in 
the development of NSCL/P in Asian populations. 

As far as we know, there has been four published meta-analyses 
regarding MTHFR polymorphisms and NSCL/P risk. According 
to Verkleij-Hagoort et al.(2007) [19] and Johnson and Little (2008) 
[20], no significant associations between MTHFR C677T and 
A1298C polymorphisms and NSCL/P were acquired. The pooled 
results of Luo Y et al. (2012) [17] indicated that maternal MTHFR 



677TT genotype was related to increased risk of having a NSCL/ 
P offspring. The more recently meta-analysis conducted by Pan Y 
et al. (2012) [18] showed that MTHFR C677T polymorphism 
contributed to elevated risk of NSCL/P among Asians. Compared 
with the previous meta-analyses, we updated this meta-analysis by 
adding newly published studies, which were not included in the 
previous meta-analyses. Finally, we achieved consistent conclu- 
sions with Luo Y et al. (2012) and Pan Y et al. (2012). 

Some limitations of this meta-analysis should be acknowledged. 
Firstly, the sources of control among the studies were different 
from each other. Some studies were population-based studies, and 
others were hospital-based studies. Secondly, our results were 
based on unadjusted OR values that lack the original data from 
the eligible studies, which could lead to relatively weak power to 
estimate the real relationship. Thirdly, our analyses were based on 
single-factor estimates, which overlooked the interactions of gene- 
gene and gene-environment in the development of NSCL/P. 
Finally, the sample size was relatively small to investigate the 
association between MTHFR C677T and A1298C polymor- 
phisms and NSCL/P risk. 

In conclusion, our meta-analysis suggested that MTHFR 
C677T polymorphism was related with an increased NSCL/P 
risk in Asian populations. In the future, large sample studies should 
be warranted to investigate the association between MTHFR 
C677T and A1298C polymorphisms and NSCL/P, and to 
examine the potential gene-gene and gene-environment interac- 
tions. 
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